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Abstract 

This  paper  concerns  the  existence  of  interned  solitary  waves  moving 
with  a  constant  speed  at  the  interface  of  a  two-layer  fluid  with  inflnite 
height.  The  fluids  are  immiscible,  inviscid,  and  incompressible  with  con¬ 
stant  but  different  densities.  Assume  that  the  height  of  the  upper  fluid 
is  infinite  and  the  depth  of  the  lower  fluid  is  finite.  It  has  been  formally 
derived  before  that  under  long- wave  assumption  the  first-order  approx¬ 
imation  of  the  interface  satisfies  the  Benjamin-Ono  equation,  which  has 
algebraic  solitary-wave  solutions.  This  paper  gives  a  rigorous  jDroof  of 
the  existence  of  solitary-wave  solutions  of  the  exact  equations  govern¬ 
ing  the  fluid  motion,  whose  first-order  approximations  are  the  algebraic 
solitary-wave  solutions  of  the  Benjamin-Ono  equation.  The  proof  relies 
on  estimates  of  integral  operators  using  Fourier  transforms  in  T^(R)- 
space  and  is  different  from  the  previous  existence  proof  of  solitary  waves 
in  a  two-layer  fluid  with  finite  depth. 
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1.  Introduction 


This  paper  deals  with  internal  waves  of  finite  amplitude  in  a  stable  two-fiuid  system 
with  infinite  height.  The  two  fluids  are  immiscible,  inviscid,  and  incompressible  with 
different  but  constant  densities.  The  lower  boundary  is  horizontal  and  no  upper  boundary 
exists.  Let  the  depth  of  the  lower  fluid  be  h  in  an  equilibrium  state  and  the  densities  of 
the  lower  and  upper  layers  of  the  flow  be  respectively,  with  p~  >  p'^  >  0.  It  is  of 

interest  to  find  an  internal  solitary  wave  moving  with  a  constant  velocity  U  at  the  interface. 
This  problem  was  first  studied  by  Benjamin  using  long- wave  theory  [4].  In  the  long- wave 
analysis,  the  lengths  of  waves  at  the  interface  are  assumed  to  be  considerably  lafger  than 
/i,  but  no  long-wave  approximation  is  possible  for  the  motion  of  the  upper  fluid  since  the 
height  of  the  upper  fluid  is  infinite.  In  reference  to  an  (a:*, y*)-coordinate  system  moving 
with  the  same  speed  U,  the  fluid  motion  is  steady.  It  is  found  in  [4]  that  to  have  long  waves 
at  the  interface,  U  must  be  near  a  critical  wave  .speed  Uq  =  ((p~  —  p'^)gh/ p~y^'^ ,  where 
g  is  the  gravitational  acceleration  constant.  Assume  that  (Uq/U)^  =  1  +  •y\€  with  e  >  0 
a  small  long- wave  pai’ameter  and  71  <  0  a  fixed  constant  and  the  long- wave  expansion  of 
the  solution  in  the  lower  fluid  is  used.  A  model  equation,  called  Benjamin-Ono  equation, 
can  be  derived  as  follows. 
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and  g*{x)  is  the  deviation  of  the  interface  from  the  equilibrium  state.  The  equation  (1) 
has  a  solitary-wave  solution 
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where  a  —  —(471/3),  p  =  —{p’^j p~'yi )  and  x  is  a  stretched  horizontal  variable  with  x  —  ex*. 

The  Benjamin-Ono  equation  was  first  derived  in  [4]  for  very  general  stratified  fluids 
with  infinite  depth  and  the  two-fluid  flow  is  one  of  these  cases.  Also  the  algebraic  solitary- 
wave  solution  (2)  was  obtained  using  Fourier  transform  technique.  At  the  same  time,  Davis 
and  Acrivos  [11]  obtained  the  equation  in  an  experimental  and  numerical  study  of  the  wave 
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phenomena.  Then  Ono  [19]  derived  a  corresponding  time-dependent  equation,  which  weis 
later  studied  using  inverse  scattering  method  [9,  12]  and  has  multi-soliton  solutions  [17]. 
Maxworthy  [18]  carried  out  experiments  again  to  show  the  existence  of  internal  solitary 
waves  and  found  that  the  waves  are  quite  stable.  Maslowe  and  Redekopp  [16]  generalized 
the  Benjamin-Ono  equation  to  nonlinear  waves  in  stratified  shear  flows  and  obtained  more 
complicated  integral  model  equations.  Then  Benjamin  [6]  considered  a  two-layer  fluid 
with  large  surface  tension  at  interface  and  derived  the  equation  (1)  with  an  extra  term 
of  third-order  derivative.  Recently  Amick  and  Toland  [2,  3]  proved  that  the  solitary-wave 
solution  (2)  and  its  translations  in  the  r-direction  are  the  only  solutions  of  (1)  which  decay 
to  zero  as  x  approaches  infinity.  Bennett  et  al.  [8]  obtained  the  stability  of  solitary  waves 
of  the  Benjamin-Ono  equation  with  time  dependence  by  a  similar  technique  used  in  the 
proof  of  stability  of  solitary  waves  for  the  KdV  equation  [5,  10].  However,  we  note  that 
the  solitary- wave  solution  (2)  is  a  solution  of  the  Benjamin-Ono  equation,  which  is  an 
approximate  equation  of  the  exact  equations  governing  the  fluid  motion.  It  is  importment 
to  know  whether  the  solitary-wave  solution  (2)  is  an  approximation  of  a  solution  of  the 
exact  equations.  The  corresponding  problem  for  two-fluid  flows  with  finite  depth  has  been 
studied  intensively  in  last  fifteen  years.  The  model  equation  for  the  interface  is  so-call 
KdV  equation,  which  also  has  solitary-wave  solutions.  It  has  been  proved  by  numerous 
hteratures  that  the  solitary-wave  solution  of  the  KdV  equation  is  an  approximation  of 
a  solitary-wave  solution  of  the  exact  equations.  The  reader  is  referred  to  an  article  by 
Benjamin  et  al.  [7]  for  a  general  treatment  of  such  problems  and  the  references  therein.  As 
pointed  out  in  [7],  the  question  whether  the  solitary- wave  solution  (2)  is  an  approximation 
of  a  solution  of  the  exact  equations  for  a  two-fluid  flow  with  infinite  depth  needs  to  be 
answered. 

This  paper  gives  a  rigorous  proof  of  the  existence  of  solitary  waves  at  the  interface 
for  a  two-layer  fluid  with  infinite  depth,  whose  first  order  approximations  are  the  solitary 
waves  given  by  (2).  The  result  can  be  stated  as  follows.  Assume  that  the  depth  of  lower 
fluid  is  h  at  infinity  and  the  height  of  the  upper  fluid  is  inflnite.  The  densities  are  p'^  and 
p~  for  the  upper  fluid  and  lower  fluid,  respectively.  If  the  wave  speed  U  on  the  interface 
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satisfies 

U  =  (|(/)“  -  p'^)gh/{p-{l  +  7ie)))  ^ 

with  7i  <  0  a  fixed  constant  and  6  >  0  a  small  parameter,  then  there  is  a  solitary-wave 
solution  for  the  exact  equations  governing  the  two- fluid  flow.  Assume  h  +  77*(x*)  is  the 
distance  between  the  rigid  bottom  and  the  interface.  Then  the  form  of  the  interface  for 
the  solitary- wave  solution  satisfies 

where  a  = -(471/3),  jU  = -(/)+//5“7i), 
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y  sup  1(1  +  x^~)(d^Y*{x)/dx^)\  <  Ke, 

and  K  is  independent  of  e  but  may  depend  on  s.  Therefore  the  solitary-wave  solution 
(2)  of  the  Benjamin-Ono  equation  (1)  is  the  first-order  approximation  of  a  solitary-wave 
solution  of  the  exact  equations. 

The  idea  of  the  proof  can  be  smnmarized  as  follows.  Since  the  interface  of  the  two 
fluids  is  part  of  the  solution,  the  domain  of  the  fluids  is  unknown.  We  use  Levi-Civita’s 
method  to  transform  the  unknown  domain  into  a  fixed  domain  using  velocity  potential 
and  stream  function  as  independent  variables.  The  main  advantage  of  Levi-Civita’s  trans¬ 
formation  is  to  put  all  the  nonlinear  terms  in  the  boundary  conditions  rather  than  in 
the  differential  equations.  It  has  been  shown  that  such  transformation  is  one  to  one  and 
onto  if  the  amplitude  of  the  solutions  is  small  [20].  Then  we  use  Fourier  transform  to  in¬ 
vert  the  differential  equations  with  nonlinear  boundary  conditions  into  integral  equations. 
Thus  finding  a  solution  of  the  differential  equations  becomes  looking  for  a  fixed  point  of 
corresponding  integral  operators.  However,  the  kernels  in  the  integral  operators  have  a 
singularity.  In  order  to  remove  the  singularity,  a  solvability  condition  must  be  satisfied. 
Next  we  introduce  an  unknown  function  and  obtain  an  integro- differential  equation  for  this 
function  from  the  solvability  condition.  By  using  some  techniques  in  operator  theory,  we 
are  able  to  transform  the  integro- differential  equation  into  an  integral  equation.  By  taking 
advantage  of  the  fact  that  the  nonlinear  terms  only  appear  at  the  boundaries  and  obtain¬ 
ing  the  estimates  of  the  integral  operators,  we  can  show  that  the  operators  are  contraction 


4 


in  certain  Banach  spaces.  The  contraction  mapping  theorem  implies  the  existence  of  the 
solution.  We  note  here  that  the  special  forms  of  the  equations  in  the  two-fluid  flow  are 
important  to  our  approach.  The  method  also  works  without  many  changes  for  a  two-fluid 
system  with  finite  height  in  upper  fluid  and  infinite  depth  in  lower  fluid,  which  has  been 
studied  by  Amick  [1]  (it  was  known  by  the  author  after  the  major  work  in  this  paper  was 
completed).  Amick  used  a  different  formulation  of  the  problem  and  a  different  approach  as 
well  to  prove  the  existence  of  the  solitary  wave  solutions.  Kirchgassner  [15]  also  considered 
surface  waves  on  a  layer  of  fluid  with  infinite  depth.  However,  we  note  that  the  ai'gument 
given  here  does  not  apply  to  the  proof  of  the  existence  of  solitary  waves  in  continuously 
stratified  fluids  with  infinite  depths,  which  were  included  in  [4],  and  further  research  is 
needed  for  such  cases. 

This  paper  is  organized  as  follows.  In  Section  2,  all  the  variables  are  nondimensional- 
ized  and  Levi-Civita’s  method  is  used  to  transform  the  governing  equations  in  an  unknown 
domain  into  equations  in  a  fixed  domain.  In  Section  3,  a  formal  asymptotic  expansion  of 
solutions  is  used  and  solitary- wave  solutions  are  formally  derived  using  a  different  method 
from  the  one  used  in  [4].  In  Section  4,  the  equations  in  Section  2  are  transformed  into  inte¬ 
gral  equations  and  a  long- wave  scaling  for  the  horizontal  variable  is  assmned.  A  solvability 
condition  is  introduced  to  remove  the  singularity  of  the  kernels  in  the  integral  operators. 
In  Section  5,  the  estimates  of  the  integral  operators  are  obtained  using  properties  of  the 
Fourier  transform  in  X^(R).  Finally,  the  existence  of  solitary- wave  solutions  is  proved  in 
Section  6. 

2.  Formulation 

In  this  section,  we  use  a  formulation  devised  by  Levi-Civita  for  treating  the  problem 
of  existence  of  waves  of  finite  amplitude  [14,  20].  The  method  is  motivated  by  the  desire 
to  reformulate  the  problem  in  terms  of  the  velocity  potential  ip*  and  stream  function  ■0* 
as  independent  variables  in  order  to  work  in  a  fixed  domain  in  stead  of  in  the  partially 
unknown  domain  in  the  physical  plane.  We  consider  an  irrotational  flow  of  two  immis¬ 
cible,  inviscid,  and  incompressible  fluids  of  different  but  constant  densities  bounded  by  a 
horizontal  rigid  bottom  and  without  upper  boundary.  There  exists  an  interface  between 


these  two  fluids  and  there  is  a  wave  at  the  interface  moving  with  a  constant  speed  U.  In 
the  reference  to  a  coordinate  system  moving  with  the  wave  at  the  same  speed,  the  flow 
can  be  regarded  as  a  steady  flow.  Let  y*  —  0  he  the  interface  when  x*  goes  to  infinity, 
y*  —  —h  be  the  rigid  bottom  with  h  >  0,  and  y*  =  rj*{x*)  be  the  interface.  Then  the 
upper  fluid  with  density  p"*"  is  in  'rf{x*)  <  y*  <  +oo  and  the  lower  fluid  with  density  p~ 
is  in  —h  <  y*  <  p*(x*),  respectively. 

Since  the  flow  is  irrotational  and  the  fluid  is  incompressible,  we  can  introduce  the  ve¬ 
locity  potential  <p*,  the  stream  function  rp*,  and  the  comi^lex  velocity  potential  X*{^*  iV*)  = 
X*{z*)  with  X*  =  +  hp*,  z*  =x*  +  iy*  in  -h  <  y*  <  r]*{x*)  and  rf{x*)  <  y-*  <  -foo. 

The  complex  velocity  w*  =  dx* /dz*  is  given  by  w*  =  u*  —  iv*,  where  u*,  v*  are  the  velocity 
components.  The  corresponding  dimensionless  quantities  will  be 

z  =  zyh,  w  =  w'/U,  p  =  p+/p-<l,  x  =  Vi+i’/>l=x7(W)  (3) 

SO  that  dx/dz  —  w.  Furthermore  we  introduce  a  nondimensional  parameter 

7  =  -  p)9h/U^  ,  (4) 

where  is  so-called  nondimensional  wave  speed.  Therefore  ipi  is  a  constant  along 

the  interface  and  for  the  sake  of  convenience  choose  i/>i  =  0  at  the  interface.  The  bottom 
corresponds  to  V’l  =  —  1  since  the  total  flow  over  a  curve  extending  from  the  bottom  to 
the  interface  is  Uh.  In  the  following  we  use  to  denote  a  quantity  /  in  the  upper- 

and  lower-layer  fluid,  respectively.  The  boundary  conditions  are  formulated  as  follows: 

V  =  —  Im  to  =  0  at  •01  =  —1,  (5) 

jy  (l/2)|to“p  —  (p/2)|t(;'^P  =  co7istant,  (6) 

y+  =  y-  at  01  =  0,  (7) 

where  the  second  condition  is  from  the  Bernoulli’s  equation  on  the  interface,  w'^,y'^  and 
w~,y~  at  01  =  0  are  denoted  the  values  of  to,y  as  0i  — »  At  the  infinity,  we  have 

1  when  |x|  — »•  -f  oo  .  (8) 
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Therefore  the  complex  potential  function  xi^)  niaps  the  physical  plane  (i.e.  the  z-plane) 
onto  the  strip  regions  —  1  <  •t/’i  <  0, 0  <  t/’i  <  +oo  in  the  x-plane.  The  mapping  is  one  to 
one  and  the  inverse  mapping  function  z{x)  exists  as  we  can  easily  see  later.  Taking  x 
independent  variable,  we  then  determine  the  analytic  function  iy(x)  in  —  1  <  V)!  <  0,0  < 
?/>i  <  +00  from  the  boundary  conditions. 

It  is  convenient  to  simplify  the  form  of  the  boundary  conditions  by  following  Levi- 
Civita’s  idea  [14,  20]  to  replace  the  dependent  variable  w  hy  6  -\-i\  through  an  equation 

w  =  exp  ^  —  i{9  +  .  (9) 

Therefore  |t<;|  =  and  6  =  arguJ,  where  |rD|  is  the  magnitude  of  the  velocity  vector 
and  0  is  the  inclination  relative  to  the  a;-axis.  Now  we  formulate  the  conditions  for  the 
determination  of  9  and  A  in  the  (fi-tpi  plane.  Condition  (5)  is  reduced  to 

^  =  0  at  ■01  =  •  (10) 


(11) 
(12) 

When  |(pi  I  — >  +00, 

^  +  z'A-^0.  (13) 

Since  9  iX  is  an  analytic  function  of  (pi  +  i0i ,  0  and  A  satisfy  the  Cauchy- Riemann 
equations 

9<Pi  =  A^i  >  ~  —  1  <  01  <  0,  0  <  01  <  +O0.  (14) 


By  differentiating  (6)  and  (7)  with  respect,  to  <^i  along  0i  =  0,  we  have 


dy  _  d\w 

1-z - h  w  Hs — 

5<pi  d^x 


dy  _  dy'^ 
d(px  9(pi 


at  01  =  0 . 


From  dz/dx  =  (i^  +  fi»)|t<i|  ^  and  the  analyticity  of  9  +  zA,  we  obtain  [dyld^px) 
and  (5A/5<pi)  =  — (5^/50i).  Thus  at  0i  =  0  (6)  and  (7)  become 

2A-  d9~  2A+  d9'^  _x-  .  - 

—--pe^^  - - je  ^  sm9  =0, 

001  atpx 

e~^  sm9~  =  sin^'*' . 
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In  this  paper  we  are  interested  in  finding  bounded  nontrivial  solutions  (^,A)  of  (14)  with 
the  boundary  conditions  (10)  to  (13)  where  A  is  even  and  6  is  odd  in 

3.  Formal  derivation  of  solitary  wave  solutions 

Although  a  formal  solution  has  been  obtained  by  Benjamin  [4]  and  many  others  in 
different  ways,  we  give  another  derivation  here,  which  not  only  gives  the  first  order  ap¬ 
proximation  but  also  leads  to  higher  order  approximations. 

To  obtain  an  approximate  solution  of  (10)  to  (14),  we  use  a  new  parameter 
to  stretch  the  horizontal  variable  Then  we  also  strech  the  vertical  variable  by  7/>  =  ei/’i 
in  the  upper  fluid,  since  the  vertical  length  is  infinite  in  the  upper  fluid  and  no  long-wave 
approximation  can  be  assumed,  and  let  the  vertical  variable  in  the  lower  fluid  be  the  same, 
i.e.  V’  =  V*!-  Assume  that  6,  A  and  7  have  the  asymptotic  expansion  of  the  following  forms, 
9+  =  e'^9+(ip^rj;)  +  e^9+((p^^)  + 

in  0  <  V’  <  +00, 

A+ =  e^Aj''(vP,?/>)  +  e^\J(<^,V’)  H - , 

-f- e'’^^(v?,t/’)  H - , 

in  —  1  <  ?/>  <  0, 

A“  =  eA;^ (v7, 7/>)  -h  e^A^ {(p,  tp) -\ - , 

7  =  To  +  eyi  .  (15) 

In  terms  of  the  equation  (14)  becomes 

Aj  =  0+,  9+  =  -Aj  in  0  <  V'  <  +c«, 

-I  <rp  <0. 

Substitution  of  the  asymptotic  forms  for  9  and  A  yields  different  order  of  the  approximations 
for  the  equations  and  boundary  conditions.  The  first-order  approximation  is 

in  0  <  V'  <  +c«, 

Ar^-  =  0’  =  ~^^<P  in  -  1  <  V’  <  0 , 

-  7o6'f  =  0,  i9j“  =9'l  at  V*  =  0, 

=  0  at  =  — 1 . 
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In  order  to  have  bounded  nontrivial  solutions  A^,  we  need 


70  =  1,  iv,  ’/»)  =  +  1),  Ai  =  -S{ip), 

d+=F  F  [S^{^)]  ,  A+=F  {ik)-^\k\e-\^^^  F[S^{cp)]  , 

where  S{ip)  is  even  and  to  be  determined,  and 

f  [/]W  =  ^  ^ 

V  /‘+°° 

F  I/l(»’)  =  /  .  (16) 

J  —oo 

are  tlie  Fourier  and  inverse- Fourier  transforms.  The  properties  of  these  transforms  can  be 
founded  in  [21,  22]. 

To  find  an  equation  for  S{(p)^  we  need  second  order  approximation  of  the  equations, 


in  0  <  Ip  <  +00, 

=  0,  0-^  =  -A^^  in  -  1  <  7/1  <  0  , 

^  -  ToSj"  =  7i«r  -  -  7oAr«r. 

9^  —  at  1/1  =  0, 

9^  —  0  at  tp  =  —1 . 


This  is  a  nonhoinogeneous  boundary  value  problem.  In  order  to  have  a  solution  (A2  ,^2)1 
the  nonhoinogeneous  terms  at  V’  =  0  must  satisfy 

(^T,sr-2Ar^+^^-7oArer)  =(i- 

By  the  first  order  approximation  and  70  =  1,  5'(<p)  must  satisfy 

■r,SM  +  3SMSvM  +  pF  -|*:|  r  [5„(¥.)|  {^)  =  0. 

Since  we  are  looking  for  solutions  decay  to  zero  at  infinity,  5((/?)  satisfies 

pF  Wf[S(y.)l  (7.)-7.S(y>)-(3/2)5"(<p)  =  0,  (17) 
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which  is  so-called  Benjamin-Ono  equation  [4,  19],  This  equation  has  a  solitaiy  wave 
solution, 


S(¥>)  = 


afi 


(18) 


where  a  =  — (471/8),//  —  with  71  <  0.  The  derivation  of  solution  (18)  from  (17) 

can  be  found  in  [4].  Therefore  the  first-order  approximation  in  (15)  can  be  obtained.  The 
first-order  approximation  of  the  wave  profile  at  the  interface  is 


V  =  <S{p)  +  0(^^)  =  -|^  +  0(6^).  (19) 

^  'T  P 

Since  a  >  0  with  71  <  0,  this  is  an  elevation  solitary  wave.  This  formal  procedure  can  be 
easily  carried  over  to  find  the  higher  order  approximations. 

In  the  following  sections,  we  shall  show  rigorously  that  there  exists  a  solitary  wave 
solution  of  (10)  to  (14)  and  (19)  is  the  first-oi’der  approximation  of  this  solution  at  the 
interface  ?/>i  =  0. 


4.  Transformations  and  Banach  spaces 

In  this  section,  we  transform  (10)  to  (14)  into  integral  equations.  Let  7  =  1  -|-7ie  with 
7i  <  0  and  e  >  0. 

First  we  consider  the  equations  in  the  domain  0  <  i/’i  <  4-00.  The  equation  (14) 
implies 

=^*,.  =-<.  =»■  (20) 

If  )(■*■,)(“  and  9~  are  small  enough,  (12)  can  be  written  as 

=  arcsin(e'^^“'^  sin^“)  =  h{(p\)  at  V’l  0-  (21) 

By  applying  the  Fourier  transform  on  (20),  we  have 

(•P  [«+])*.♦.- 1*^1"  f  [«■"!  =  0, 

A  A  ^ 

where  F  [/]  is  defined  in  (16).  If  F  [0+]  is  finite  as  ipi  +00, 

F[e+](k,^h)  =  e-\’^^^^ 
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which  gives 


F  arcsin(e^'^-^  sin^“)  .  (22) 

L  L  V'i=oJJ 

By  (20)  again,  we  have 

F  arcsin(e^^~'^  sin0~)  .  (23) 

L  ’/'i=oJJ 

Note  that  (22)  and  (23)  hold  in  0  <  i/^i  <  +oo. 

Next  we  study  the  equations  (10)  to  (14)  in  the  domain  —1  <  V’l  <0.  The  equations 
are 


=  0 

in  —  1  <  x/^i  <  0 

aV’i 

=  7ie^  +  ^  ^  ' 

Ofi 

+  (1  +  7ie)(e-^^  .sin(^“)  ■ 

-r 

)  at  '(/’i  =  0, 

o 

II 

1 

at  tpi  =  — 1 . 

(24) 

(25) 

(26) 


From  (22),  we  can  rewrite  (25)  as 

d6~  -  r  r-  T 

—  -e-+pF  |Ar|F[rJ  =7ie^- 

+  (1  +  7ie)(e“’^^"  sin(r )  -  r)  +  p(l  -  F  \k\  F  [^"1 

+  ^F  \k\  F  —  axcsin(e^^“'^  sin^“) 

=/i(r,A-,A+).  (27) 


By  using  the  Fourier  transform,  (24),  (26)  and  (27),  we  obtain 


0  {}piA\)  =F 
A“(v?i,V’i)  =F 


sinh(|A:|(t/>i  +  1))  F  [/i] 

1^1  cosh  |/:|  —  sinh  |A:|  +  p\k\  sinh  |A:|  ’ 

|fc|c.osh(|A;[(^i  +  1))  F  [/i] 
zfc(|A:|  cosh  |fc|  —  sinh  |A:|  +  p\k\  sinh  |^|) 


(28) 

(29) 


Thus  the  equations  (10)  to  (14)  are  transformed  into  integral  equations  (22),  (23),  (28) 
and  (29),  where  f\  is  defined  in  (27). 
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A  crucial  step  in  our  procedure  is  that  we  make  the  same  change  of  independent 
variables 

(f  =  eipu  ’Z'  =  ’Ai  ,  (30) 


which  we  made  in  Section  3.  In  accordance  with  the  development  of  Section  3,  we  put 


0±  =  e-^e^,  A±  =  e-^A*. 


±  _  .-1  \± 


(31) 


Then  by  changing  the  integration  variables  in  (22),  (23),  (28)  and  (29),  these  equations 
become:  In  0  <  i/’  <  +oo. 


^+((^,^)  =F 


A+((^,V>)  =F 


jp 


€  ^  arcsin  ^sin(e^^  ) 


^—0 


arcsin  ^  sin(e^^  )  )  ;  (33) 


In  — 1  <  V’  <  0, 


9  =F 


X-M)  =F 


sinh(e|A:|(?/>  +  1))  F  [/2] 


pe|A:|  sinhe|A:|  +  e\k\  cosh e\k[  —  sinh e|ifc| 

_ |A:|  cosh(6|fc|(i/i  +  1))  F  [e/2] _ 

ik(pe\k\  sinhe|A:|  +  t\k\  coshe|A:|  —  sinhe|A:|) 


where 


(32) 


(34) 

(35) 


M0-,X-,X+)  =  6-2 Me^9-,eX-,  eX+) 

=  [l\^9~  4-  (1  +7ie)(e“^e“^''^  sin(e^6'~)  -  ^~)] 
+  F 

+  p{l 

=  J  +  II  + JJJ. 


e\k\  F  9  —  e  ^  arcsin(e' 

6|A:1  F  [r] 


iA+-eA' 


sin(e20-)) 


Notice  that  when  e  or  A;  is  small,  the  denominators  in  the  inverse-Fourier  transforms  of  (34) 
and  (35)  have  a  zero  point  of  order  e^k'^ .  Therefore  we  have  to  single  out  this  singularity 
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in  the  integrals.  Write 


sinh(e|fc|(^  +  1))  _  'ip  +  1 

p€\k\  sinhe|A:|  +  e|A:|  coshe|fc|  —  sinhe|A;|  pe|Ar| 

('i/>  4-  l)(e|fc|  cosh(e|^|)  —  sinh(e|fc|)) 
p€\k\(p€\k  \  sinhe|A;|  +  e|A:|  coslie|fc|  —  sinhelA;|) 

^  sinh(e|fc|('0  +  1))  —  (V’ +  1)  sinh(e|fc|) 
pe\k\  sinhe|A:l  +  €|fe|  coshe|A:|  —  sinhe|fc| 

ki{k){tp  +  1)  +  k2{k){tp  +  1)  +  kz{k,%p),  (36) 


|A:|  cos^el^Ki/j  +  1))  1 

ik{pe\k\  sinhe|^|  +  e\k\  coshe|A;|  —  sinhe|fe|)  ipe'^k\k\ 
(e|fc|  cosh(e|^|)  —  sinh(e|A:|)) 
ip€'^k\k\{pe\k\  sinhe|A:|  +  e|fc|  coshe|fc|  —  sinhg|A:|) 
(e|A;|  cosli(e|^|(i/>  +  1))  —  sinh(e|A;|)) 
ike{pe\k\  sinhe|A;|  +  e\k\  coshe|A:|  —  sinhe|A:|) 

{1  / iek)  ^ki{k)  +  k2ik)  +  k3^{k,rp)^. 

Therefore  (34)  and  (35)  become 


{l/p€\k\)  F  [f2] 


+  F 


h{k)  F  [f2] 


+  F 


k,(k,i^)  F  [f2] 


A-(vp,V>)  =F 


{l/ip€k\k\)  F  [/s] 


+  F 


{l/ik)k2{k)  F  [f2] 


+  F 


{l/ik)k,4k,i>)F[f2] 


(37) 


(38) 


(39) 


Note  that  ks(k,'ip)  is  finite  as  ek  0.  Therefore  we  need  to  remove  first  two  terms  in  the 
right  hand  sides  of  (38)  and  (39).  Let 


A  =  -r]{(p)  +  C  and  then  B  =  O’  +  ^  ■ 


(40) 


If  we  assume 


—T]  =F 


{l/ipek\k\)  F  [/a] 


+ 


V 

F 


{llik)k2{k)  F  [f2] 


(41) 
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which  is  called  solvability  condition,  then  (38)  and  (39)  are  changed  into 


C=F 


h{k,i^)F[f2] 


{l/ik)k3^{k,7jj)  F  [f2 


^=0 


By  the  definition,  ,A  jA"*")  can  be  written  as 

f2(0-,X-,  A+)  =  f2iv<p  +  +  C,  A+) 

=  71  er  -  zex-e-  +  e^n{e-,x-,  a+) 

=  7i  +  3er/T;<^  +  7i  -  HC^^p  “  <^7  +  CO 

+  e27^(^",A",A+) 

=  7ie7<^  +  C,  C.  A"*”) . 


Thus  by  taking  the  Fourier  transform  on  both  sides  of  (41),  we  can  obtain 


-pe\k\ki  F  [7/]  =  (1  +  pe\k\k2{k))  F  [/a] , 


which  implies 


F 


p\h\(l  + /x\k%(k))  '  F[p] 


=  (l/e)/2(r,A-,A+) 


J  (p 


iiVp  +  +  Mv,  C,  C,  A+). 


If  we  let  e  =  0  and  /2  =  0,  (45)  becomes  (17)  and  it  has  a  solution 


5((^)  = 

where  a  =  ’-'(471 /3), =  “(p/Ti)*  Assume 


(^2  -I-  i^i2  ' 


7  =  S{ip)  +  uj{p). 


Then  substitution  of  (47)  into  (45)  yields 


plF 


\k\{l  +  pe\k\k2{k))  ^  F  [u 


-  -  S{Suj)^ 


—  (3/2)(u)  )^  +  p  [F 


\k\^pek2{k){l  +  pe\k\k2(k))-^  F  [5(¥p)] 


+  f2i 


By  applying  tlie  integration  once,  we  have 


pF  \k\{l pe\k\k2{k))  ^  F  [to]  —'Jiuj  —  SSuj 

V  f  A  1  r+oo 

=  {3/2y +  pF  \k\^pek2{k)il  +  pe\k\k2{k))-^  F[S{<p)]  + 

=  (48) 

In  tlie  following,  we  shall  show  that  (32),  (33),  (42),  (43)  and  (48)  have  a  solution 
(cu,  1^,  C,  A"*")  in  some  Banach  space  for  small  e  >  0,  where  u;,C,A'''  are  even  and  i^, 

j 

axe  odd  in  t 

Before  proving  the  existence,  we  define  some  Banach  spaces  for  later  use.  Let  jy^(R) 
be  the  usual  Sobolev  spaces  on  R  =  (— oo,  +oo)  with  5  >  0  [22].  Define 

W^  =  {f{^)eH^{K)\cpf{^)eH%K)  and  ||/|k^  =  ||(1  + 

=  ifiVii’)  ^  L^(R  X  [-1,0])  for  any  ■)/>  G  [-1, 0], /((p, V’)  G  W^,  and 
ll/ll S-*  =  sup  ||/((^,t/j)|k..  < +oo}  . 

We  use  IF/,  Bl,Hg  or  IF/,  Bl,Hg  to  denote  the  corresponding  spaces  with  even  functions 

A 

or  odd  functions  in  (p.  Note  that  f{(p)  G  if  and  only  if  F  [f]{k)  is  differentiable  in  k 
and 

(1  +  P)^/2  F  [/](^)  G  I'(R),  (1  +  F  [/](A;)  G  ^'(R) . 

Also  by  using  the  Sobolev  embedding  theorem,  we  have 
Lemma  1.  If  u,?;  G  H^(R)  for  s  >  1,  then  uv  G  H^(R)  and 

where  AT  is  a  fixed  constant. 

5.  Estimates  of  the  operators 

First  let  us  consider  (32)  and  (33).  Rewrite  these  equations  into 

e+ip,xl;)=F  E  [5,^]  +F  e-^l"l’^E[u»^  +  e(<^,0)] 
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(49) 


+  F 


^-e\k\rk  p 


■  ^  arcsin(e*^^'*'  ^  ^  sine^6  )  —  6' 


l/;=:0 


A+(9,«  =f 

(ii/ble/fc)e-^l*'^  F  [5,^] 

V 

+  F 

(i|A;|e/A;)e-*l'=l^  Ak+e(v^,0)] 

V 

+  F 

(ee|fc|/A:)e-^l^l^  F 

arcsin(e^^^^~^  ^  sin  e^^"")  — 

V'=Oj  J 


def 


Ml  (5)  +  M2(a;,  0  +  MjCo;,  C,  A+), 


(50) 


where  Lz,Mz  are  of  order  e  for  small  e  >  0.  By  using  the  properties  of  Fourier  transforms, 
Lemma  1,  and  the  Plancherel  equality,  we  can  have  that  if  |||(v^,0)||w«  +  |iC(‘P) 0)11  + 

||a;<^||i4/3  +  ||A+((^,0)||vv>  <  K  for  .s  >  2,  then  for  0  <  V*  <  +oo, 

||L2(td,0IU,'  +  ||M2(a;,e)||//|  <  A^KIk^  +  ||e(¥^, 0)|k^ ) 


and 


||A3(w,(f,C,  A'*')||//»  +  ||M3(cd,^, C,  A^)||/f|  <  Ke. 


Since  I  A:  I  is  piecewise  differentiable  and  are  odd, 

r-\-co 

(pM2{oj,^)  =  2ip  cos(Ar£p)(ie)e“*l*l^  F  [u^  +  ^(y’jO)]  dk 
Jo 

=  -2  /  sin(k^)(u)  (^e-'*'*  F  K  +  «^.0)1  j  dk 

/  pI  a  \ 

e-“*’(sgn(*;))e  (^e”*'*'*  F  +  ?(v,  0)] j  dk 
where  sgn(A:)  =  1  for  k  >  0  and  sgn(fc)  =  —1  for  A;  <  0.  Thus  for  ^  bounded, 

\\<pL2{oj,^)\\h‘  +  \\^M2{uj,0\\h’ 

<  A'(||(l  +  ky/‘^  (d  F  k  +  e(<^>  0)]^  /dk 


L2(R) 

+  I  (1  +  k’^y/hgXiik)  [d  (^(e-'l"^!^  F  K  +  e(‘^, 0)]  )  /dk 
<A-(|Kk. +||e(<^,0)|k,). 

By  using  the  same  argument,  we  have 

?,  c,  A+)||„,  +  \\vM,{u,,  J,  C,  A+)||„.  <  Ke 


Z,2(R)/ 
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Therefore  we  obtain 

Theorem  1.  If  |iw<^||iv/  +  0)111^'  +  ||C(9^,  0)||h^'  +  <  K  with  5  >  2, 

then  for  0  <  "0  <  A'l  <  +oo 

||-Z^2(w,0l|vv-'  +  \W2{i^,0\\w-;  <  K{\\uj^\\w‘  +  ilC(<^,0)||w»)> 

Next  we  consider  (42)  and  (43).  From  the  forms  of  kz{k^'^)  and  kz^(k,t/)),  it  is  easy 
to  see  that  for  any  small  e  >  0,  A;  G  (— oo,  +oo)  and  — 1  <  0  <  0, 

\k3(k,'i/;)\  +  \kz^{k,^/;)\  +  \ekkz(k,il})\  +  |(eA;)“U.3^(^, 7/))|  <  Kz  <  +00,' 

where  K2  is  a  fixed  constant.  From  the  definition  of  /2  in  (34)  and  (35),  we  also  have 


kzik,i>)F[f2] 


=F 


kz{k,i>)F[I] 


+  F 


kz{k,t)F[II] 


+  F 


kz{k,xl^)F[III] 


Assume  that  at  0  =  0,  ||A'*‘||vy>  +  ||A  II IV;' +11^  llw';' +II^‘*'II»V;'  ^  K  for  5  >  2.  By  Lemma 
1,  the  forms  of  J,  II  and  III,  and  the  Plancherel  equality,  we  have  that  for  —  1  <  0  <  0 


h{h,i,)FlI\ 


Hi 


<K\\I\\H.<Kt. 


Also 


kz{k,i>)F[II] 


<  K 


F 


L2(R) 
A  r 


|A;|e  F  6  —  e  ^  arcsin(e*^^'''  ^  ^sine^0  ) 


L2(R) 


<Ke^{\\0  ||]/i  +  IjA  llm  +  IIA'^ll'ii) , 


and  for  s  >  1 


F 


dip 


kz{k,i^)  F  [JI] 


< 


<  K\\k^-\-^  F  [II]Hi.(R)  <  F||(l/e)(d^-^II/d<^)||i.(i,) 


L2(R) 
2 


\\k^kz{k,i>)F[II]\\l,^^) 


^5-1  V 

F 


\k\  F  0  —  e  ^  arc.sin(e'^^''^  ^  ^sine^0  ) 


<  Ke^{\\$~\\'fjs  +  ||A  \\\>  +  ||A'''||^i). 


2 

L2(R) 
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Tims  II  F  [A;3(fc,V’)  F  <  A'e.  By  using  similar  estimates  for  ///,  we  can  have 

II  F  [k3{k,‘tp)  F  [JJJ]]||//j  <  Ke.  Therefore  for  — 1  <  ■0  <  0 


F  \k3(k,^l^)F[h]  <  Ae. 

Hi 


Also  note  that  | fc |  is  piecewise  differentiable  and  the  derivatives  of  kz ,  /cs ,  ekkz  and  ( 1  / efe ) ^3 ^ 
with  respect  to  k  are  bounded.  Therefore  using  similar  arguments,  we  can  have 


^F  k3{k,,l>)F[h]  <K  F  {dk3{k,i^)ldk)F[h] 

L  J  L  J 

+  F  k3{Ki^)F[Ff2]  <Ke. 


Thus  F  [kzik,^;)  F  [f2]]  G  and 


F  kzik,-,^)  hfo.]  1^^,.  <Ae, 


By  a  similar  reason, 


F  (l/ik)k,^{k,i,)  F  [h]  <||(lAfce)fe<.(i,</>)er[/2l||iw 

<A1/2|U=(ii,  <ire. 


and  for  5  >  1, 

1^  F  ^{l/ik)k,^,(k,i,)  F  [m]  V)  F  [/2I  |U.(R| 

<md’-'h/dv-')\imK)<K,. 

V  A 

Thus  II  F  [{l/ik)k3^{k,ip)  F  [/2]]||/f''  <  Ae.  By  applying  the  similar  argument,  we  can 

V  A 

have  \\kp  F  [{l/ik)k3.^[k^'^)  F  [/2]]||//'*  <  Ke  since  |A:|  is  piecewise  differentiable.  Therefore 
F  [il/ik)k3,i;ik,xk)  F  [/2]]  e  and 

<  Ke. 

J 

We  summarize  the  above  results  as  follows. 


{l/ik)k3.^(k,il))  F  [f2 
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Theorem  2.  If  ^-  =  (5  +  +  i)  +  (  e  X-  =  -{S  +  to)  +  (  £  X+(cp,  0)  e  W, 

and  ||tu||j^,+i  +  ||i^||b'  +  IKIIbi  +  ||A+(</?, 0)|| vK'  <  K  with  3  >  2,  then 


h{k,i^)F[f2] 


B’ 


<  Ke , 


F 


{l/ik)h^{k,^l;)  F  [f2] 


Bi 


<  Ke , 


where  K  is  a  constant  independent  of  e  and  i/>. 
Finally  we  study  the  equation  (48).  Let 


N{u)  =  pF 


\k\{l  +  pe\k\k2(k))  ^  F  [u 


—  ■jicu  —  3Suj. 


Then  (48)  becomes 


iV(u,)  =  9,(u,,e,(.A+). 

We  first  study  the  general  nonhomogeneous  equation 

N{ui)  =  g{tp),  a>^(0)  =  0, 


(51) 


(52) 


(53) 


where  u>{<^)  and  are  even  functions  in  <p.  By  applying  the  Fourier  transform  on  both 

A 

sides  of  (53)  and  using  F  [5'((^)]  =  {ap,/2)  exp{— p\k\),  we  have 


p\k\{l  +  pe\k\k2{k))-^  F  H  (k)  -  71  F  [cu]  {k) 

/  +  00  ^  ^ 

-oo 

A 

where  a  and  p  are  defined  in  (46).  Thus  F  [aj](A:)  must  sati.sfy 


(54) 


F  H(Jt)- 


3ap 


2{p\k\{l  +  pe\k\k2{k))-'^  -7i) 


A 

-oo 


F  M(^) 


=  go{k)^ 


p\k\{l  +  pe\k\k2ik))  ^  -71) 

Note  that  by  the  definition  of  k2{k), 

€\k\  coshe|A:|  —  sinhe|fc| 


(55) 


(1  +  pe\k\k2{k))  =  1 


pe\k\ sinhe|A:|  +  ^\k\  coshe|A:|  —  sinhe|A:| 


-1 


^  ^  €\k  \  coshe|fc|  —  sinhe|fc|  ^  ^ 
p|fc|esinhe|^)|  ~ 


(56) 
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Denote  the  integral  operator  in  (55)  by 


TMk))  = 


Zaii 


r+oo 


2{p\k\{l  +  pe\k\k2{k))  ^  -  71)  y_c 


u{k')dk'  (57) 


def 


and  r(u)  =  Tq{u).  Then  (55)  becomes 


F[u]-%{F[u])^go{k). 


(58) 


Now  consider  an  integral  equation 


u 


T(u)  =  /. 


(59) 


where  /  G  r^(R)  and  «,/  are  even  in  (p.  Let 

/-|-oo 

t->^\'^-'^'\u{k')dk'. 

-00 

By  (59),  W{k)  satisfies 

Wkk-IJ^WF  W  =  f{k).  ■  (60) 

p\k\-ji 

Since  u{k)  and  f{k)  are  even,  we  consider  equation  (60)  with  W^i;(0)  =  0.  It  is  easy  to  find 
a  solution  ui(k)  of  the  corresponding  homogeneous  equations  (60)  for  A;  >  0, 

ui{k)  =  k{k  +  . 


Then  using  reduction  of  order,  we  can  obtain  another  linearly  independent  solution  U2{k) 
with  ^2(0)  =  0  ~  ^1^2  =  1-  Also  a.s  k  +00,  U2{k)  and  its  derivatives  have 

order  of  0(k~‘^e^^).  Thus  the  solution  of  (60)  is 

pk  ^4*00 

W(k)  =  ui{k)  U2{.s)f(s)ds+U2{k)  ui(s)f{s)ds 

Jo  J  k 

/*+oo 

N‘{f{k))'^^  I  ker{k,.s)f{s)ds  for  A;  >  0 ,  (61) 

Jo 

and  evenly  extended  to  A;  <  0. 

Lemma  2.  If  f(k)  G  Lg(R),  then  J\f{f)  G  H^(R)  and 

.  <  mhim . 
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Proof:  Since  ^{f)  is  even, 

r+oo 

lW)lll2(R)  =4|W)||^2(r,+)  <4  ker{k,s)f{s)ds 

r-jroo  /  /*+oo  \  /  r+oo  \ 

<4y  ij  I ker(A:, .s)|rfsj  yj  \keT{k,s)\f‘^{s)dsjdk 


2 

LHK+} 


A  4-00 


<4/^  sup  /  |ker(A:,.s)|rfs^  r  sup  /  | ker(fc, |  ||/||i2(R.+  ) 

Vo<fc<+ooyo  /  \0<s<+ooyo  / 


But 


sup  [  |ker(A:,s)|c?.s  <  sup  f  /  |ni(A:)u2(s)|(is 

0<k<-\-oojQ  0<k<-\-co^Jo 

+  \u2{k)ui{s)\d.s)  <  K  sup  (  [  (k^ +  l)e-^'^{s^ +  l)-^e>^^ds 

A+OO  . 

+  J  {k^  +  +  l)e~>^^d.sj  <  K,  <  +00  . 


The  same  estimate  holds  for  supo<s<4.oo  I  ker(A:,  s)\dk.  Thus  W^ffWiKK)  <  KWJWlkk)- 
By  differentiating  (61)  once  to  have 

1  pk  A-i-00 

^Afifik))=u[(k)J  U2{s)f{s)d.s +U2{k)  ui{$)f{s)ds, 

and  using  the  same  argument  for  J\f{f{k)),  we  obtain  \\d{M{f)) / dk\^i2{^-^  <  A'’||/||/,2(i{.). 
Since  A''(/)  satisfies  (60),  lld^(A/'(/))/dA:^||^2(R)  <  A'||/||/,2(r.).  Therefore  \\j'^ /Wh^CR.)  < 
Al|/||Li(R)- 

Thus  the  solution  of  (60)  is  W{k)  =  for  /  G  I'g(R).  By  (57)  and  the  definitions 
of  T  and  W{k)^  we  can  write  (59)  into 

Zafj? 


u{k)  + 


or 


u{k)  =  m- 


pifci  -71 


J^if)  =  /, 


-A/'(/)  • 


(62) 


p\k\-ji' 

From  (62),  Lemma  2,  and  the  fact  that  the  first  order  derivative  of  |A:|  is  piecewise  con- 
timxous,  we  have  that  (59)  has  a  unique  solution  u  G  ■Lg(R)  with  ||'u||l2(i{,)  =  ||(/  — 

t)-VI1l3(r)  =  \\Wkk  -  f^^wWmR.)  <  a||/IIl^(r)  and 

11(7  -  T)“V||hi(R)  =  ||w||//HR)  <  A"||/||hi(r), 

\\k{i  -  t)“VII//hr-)  -  P'^IIhhr)  -  -^(P(^//‘^^)IU2(R)  +  ||/||//i(R))- 
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Next  let  us  consider 


u  -  Te{u)  =  /  (64) 

for  e  >  0  small,  where  /  and  u  are  even  in  with  /  €  We  shall  show  that  I  —  T^ 

is  invertible  in  ie(R)  for  /  6  Ll(R).  Write  (64)  as 


u 


-r(«)  +  (T-T, )(«)  =  /. 


(65) 


Now  we  estimate  the  operator  norm  of  T  —  T^.  If  u  G  iie(R),  then  by  (56) 

3a//p|A:|(e|A:|  coshe|fc|  —  sinh  e|A:|) 


def 


Z(u)  =  iT,-T)(u)  = 


2pe|A:|sinhe|^|(/p|A:|  -  ji){p\k\(l  +  pe\k\k2{k))-'^  -  71) 

/  +  00 

e-^^l’^-'^'\u{k')dk' . 

-00 


Therefore 


(fcp(e|A;|  cosh  e\k\  —  sinhej^l)' 


/+00 

-00  (e|^|  sinhe|A;|)2(plA;|  -  7i)2(p|A;|(l  +  pe\k\k2{k))-^  -  71)2 

<  F(e|A:|  cosh  e|A;|- sinh  e|A:|)2 


(|A:|  +  l)^(e|A;|  sinhe|A;|)2 


■dk 


However 


e\k\  coshe|A;|  —  sinhe|A:| 


<  Kz  <  +00  for  all  e  >  0,  ^  G  (“OOj  +00)  • 


ej^l  sinhe|^| 

By  the  dominated  convergence  theorem,  /(e)  — »•  0  and  ||7^||  <  KI{e)  0  when  e  — >  0.  For 
small  e  >  0  the  equation  of  u  —  Te(u)  =  /  can  be  transformed  into 


u 


{I-T)-\T,-T)u  =  {I-T)-^f. 


Since  (I  —  T)~^  is  bounded  from  Tg(R)  to  Ll(R),  the  operator  norm  of  {I  —  T)~^  (T^  —  T) 
goes  to  zero  as  e  — )•  0,  which  implies  I  -  {I  -  T)~'^{Te  -  T)  is  invertible  and  the  solution 
of  a  —  Tg(«)  =  /  satisfies 

II^^IUkr)  =  IK-^- Te)“V|U|(R) 


< 


(l-(I-Tr\T.-T)y\l-T)-'f  <A1|/||,|„„,  (66) 
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if  6  >  0  is  small.  By  (64)  and  the  definition  of  in  (57),  we  also  have  that 


||^w||i,2(R.)  <  ||A:Tf(u)||i2(R)  +  P/1|l2(R)  <  -?^(P/1U2(R)  +  ||/|U2(R))  •  (67) 


For  f(k)  G  iir^(R),  by  taJking  the  derivative  of  (64),  we  obtain 


d  ...  id 

—  u(fc)  +  k—u{k) 

dk  l2(r)  dk  /-^(R) 


id/(i) 


L2(R)  II  dk  IIl2(R) 


+  IK*;"  +  1)’'V(*;)IU=(II) 


If  (1  +  k'^y/'^  f{k)  and  (1  +  k'^Y (df  f dk)  are  in  I(^(R),  then  by  using  (64)  repeatedly  and 
noting  that  the  integral  operator 

/+0O 

-OO 

maps  I'^(R)  to  I'^(R)  continuously  for  any  5  >  0,  we  obtain 


||(i  +  P)-/=«(mii.(R,+  {1  +  k^y/^ -u{k) 

<  K  (iKi  +  ky/^f{k)\\LH«.)  +  ||(i  + 


In  order  to  have  (69)  we  have  used  the  fact  that  u  =  Te(u)  +  /  implies 


l|^^^llz,2(R)  <  ll^7'e(u)P2(R)  +  ||A;/||l2(r)  <  if(p||z,2(R)  +  P/||l2(r))  . 

Now  we  can  consider  (58).  From  above  derivation,  for  e  >  0  small  enough,  (58)  has  a 

A 

solution  F  [a;]  such  that 

11(1  F  +  ||(i  + 

if  the  right  hand  side  is  finite.  Then  we  use  the  Plancherel  equality,  the  properties  of  the 
Fourier  transform,  and  the  definition  of  go{k)  in  (55)  to  obtain  that  for  g{kp)  G  with 

s  >  1,  the  equation  (53)  has  a  solution  lo  N~^{g)  €  IF/  satisfying 
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We  summarize  it  as  follows. 

Theorem  3.  If  N(uj)  is  defined  in  (51)  and  g{<f)  6  for  3  >  0,  then  for  small  e  >  0  the 
equation  iV(w)  =  g  with  a;<^(0)  =  0  has  a  solution  u  =  N~'^{g)  in  satisfying 

<K\\g\\w^. 


6.  Existence  Proof 

We  consider  (42),  (43),  (49),  (50)  and  (52).  Let  ^ 

p{(p)  p-iip)=  (^■^(v7,0)  -  Li(6'))e“^/^, 

P2{^)  =  (a+((,?,0)  -  Mi(5)^e~^/^ 

5i(¥5,'0)  =  <^((/?,  q2(p,il^)  =  (70) 


Then  we  rewrite  (42),  (43),  (49),  (50)  and  (52)  by 


=  •P(p,Pi,P2,5i,?2),  (71) 

pi(v?)  =  (X2(pe^/^,(Z1£^/^)  +L3(pe^^^,e^/^<Zi,e^''^g2,e^^V(<p)  +  -^i(-S'))) 

=  A(p,Pi,P2,5i,?2),  (72) 

P2{p)  =  {m2 (pe^ ,  5i «^ ^^ )  +  -^3  3i ,  52 ,  e^/^p2 (<p)  +  -^1  (-S')) ) 

=  A(p,Pi,P2,(ifi,<?2),  (73) 

gi(<^,./;)  =  e-'/2E  A:3(A;,V’)E  [/2((5  +  e'/^p)^(V’  +  l)  +  ^le'/^ 

-  (5  +  e^/^p)  +  ^^2e^'^^Ml(5)  +  e’/^P2) 

=  Qi(p,Pi,P2,<Zi,(Z2), 

?2(<p,^)  =  ^ 


(74) 


A  r 


{l/ik)h^ik,i;)  F  [/2((5  +  ei/"p)^(^/>  +  l)  +  ?lf'/^ 
-  (5+  e^/^p)  +  q2e^^‘^,Mi{S)  +  e^/^p2)] 


=  Q2(P,P1,P2,?1,<Z2)- 


(75) 
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We  axe  looking  for  fixed  points  of  mai^ping  (-P, -Pi ,  A,  Qi )  Q2)  in  some  Banach  spaces. 
Define  a  clo,sed  convex  set  in  a  Banach  space  x  x  W/  x  x  as  follows: 

S  ={x  =  fep,,P2,5i,52)  €  w,*+'  X  w;  X  w;  X  S*  X 

=  l|p|liv,'+'  +  IIpiIIiv;  +  ||P2||iv.-  +  Ikills;  -t-  IIkIIb;  <  <  4-oo}, 

where  61  >  0  is  a  small  fixed  number  with  s  >  2.  Assume  that  X  =  (p,  Pi ,  P2  >  <Zi  >  92 )  €  S. 


By  Theorem  2,  we  obtain 

||Qi(p,Pi,P2,5i,g2)||B2  +  ||<32(p,Pi,P2,<?i,52)|1b|  <  ^  (-70) 

In  order  to  use  Theorems  1  and  3,  we  need  to  change  (71)  to  (73)  into  the  following  forms, 

P(<P)  =^(p,Pi>P2,QuQ2)  =  P{p,Pi,P2,qi,q2),  (77) 

Pi('p)  =-Pi(^",Pi,P2,Qi,Q2)  =-Pi(p,Pi,P2,<?i,<?2),  (78) 

P2(<p)  =P2(-P,Pi,P2,Qi,<52)  =  P2(p,Pi,P2,?i,32).  (79) 

From  Theorem  3,  we  need  to  show  that 

gi{e^^^p,e^^‘^Q-i,e^^'^Q2,€^^^P2  +  Mi{S))  <Ke.  (80) 

By  (48)  and  noting  that 

|7i|  ^  F  \k\pe\k\h{k)(l  +  pe\k\h{k))-^  F  [P(<p)] 

<Ar|  F  [\k\p€\k\k2{k){l  +  p€\k\k2{k))~'^  exp(-/i|A:|)]| . 


Since  decays  exponentially  as  fc  +oo  and  k2{k){l  +  pe\k\k2{k))~^  is  uniformly 

bounded,  l|7i||pv»  <  Ke  for  any  s  >  0.  Therefore  the  W/-norm  of  first  two  terms  in  pi 
of  (48)  are  less  than  Ke.  Next  we  need  to  show  /^°°  f2{S  +  w,  C)  A'^)dcp  G  W/.  By  the 
definition  of  f2  in  (44), 

MS  +  u,  e,  C,  A+)  =  71^  +  [-3(C(5  +  cu)^  -  ^{S  +  u)  +  ce) 

+  eF{{S  +  w)(p  +  —(5  +  w)  +  C)  A*^)] 

"=  7ie  +  /i'i(u;,e,C)-  (81) 
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Now  we  need 


Lemma  3.  If  ui{(f)  G  W^^U2{^)  G  for  s  >  1,  then  ui{(p)u2{(f)d(p  G  and 


^+oo 

/  Ui{ip)u2{ip)d^  <  A"||ui((^)||w'.||u2(<(5)||w-*  . 


Proof:  Let  v{(p)  =  ui((p)u2((p)d(p.  Then  {dv/d^)  =  uiU2  G  and  by  Lemma  1 

\\{dv/dcp)\\w^  <  K\\ui{ip)\\w’\\u2W)\\w^  ■ 


Thus  we  only  need  to  show  the  estimates  for  ||v||iyo.  Since  v  is  even, 

/>+oo  ,  /•+00 


|u| 


Lim 


/•-t-oo  2 

I  ui{t)u2{t)dt  dip 
J  (p 

/»~hoo  z'+CJO 

<K  ((^^+1)“^  /  Ui{t)u2(t){l  +  t'^)dt 

Jo  J  (p 

<  A||(l  +  , 


dip 


and 


/+0O  . 

+  l)“^j  /  Ui{t)u2{t){l  +  t^)dt 

<  A'lKl  +  ip'^y^'^Ui  ||£2(R)|1(1  +  ‘y^^)^^^’«2||£,2(R)  . 


dip 


Thus  ||v||vvo  <  A'IIuiIIvk*  ||^i2||w/  and  the  proof  is  completed. 

From  (81),  it  is  easy  to  check  that  every  term  in  II\  is  either  a  product  of  ui  and 
U2  with  ui  G  kF/“^,n2  G  Wg~^  or  a  term  like  F  [|A:|  F  [f{p)]]  with  ||/||iy;  <  A'e.  The 
integral  of  the  latter  is  F  [(|fe|/i^)  F  [/((^)]]  €  W^.  Thus 

r+oo 


/ 

J  up 


III  dip 


w; 


<  K(e+  \\C\\w’  +  ll^llw^O- 


But  by  the  definition  of  P  in  (77), 


/2(5  +  a;,e'/2Qi,e^/2Q2,A+)  =  7iQiU=o+//i(u;,6'/2Qi,e'/"Q2,A+) 
=  IIiiS  +  u,e^/^Qi,e^/^Q2,\-^). 
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By  (76),  (82)  and  (83),  T  maps  S  into  itself  if  e  >  0  is  small.  Also  by  using  the  same  proof 
as  (76),  (82)  and  (83)  from  Theorems  1,  2  and  3,  we  have 

Theorem  4.  It  JTd)  =  (!-<*>, and  <;!">,??')  are 

in  S,  then 

|||T(X('))  -  T(X(2))|||  <  Ke^n\ix(^)  -X(2>|||. 

We  choose  e  so  small  such  that  the  constant  in  Theorem  4  is  also  less  than  1/2. 

Then  T  is  a  contraction  in  S.  By  the  contraction  mai^ping  theorem,  we  can  have  a  unique 
fixed  point  of  T  in  <S.  Therefore  we  obtain  a  solution  (p,pi ,P2,  (Zi,  32)  in  S  for  (71)  to  (75). 
However  from  (73),  (79)  and  (50),  it  is  easy  to  see  that  if  |inj<^||w>  +  ||^(v^,  0)||iy»  <  Ke 
and  ||A+||w/  <  K,  then  \\P2{p-,PiiP2,(Ii,q2)\\w‘  <  Therefore  we  have  the  following 

Theorem. 


Theorem  5.  For  e  >  0  small  enough,  the  integral  equations  (32)  to  (35)  have  a  solution 
e+{^,i;)=F  F  [S^]  +d+i^,i^), 

X+{ip,i;)=F  {i\k\e/k)e-^\>^\^  F[S^]  FX+M), 

=  (t/i  +  l)5<p  +  6>“(v?,'0), 

X~{ip,  -tp)  =  -S  +  X^ (9?,  Ip) , 
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where 


S(cp)  = 


aj.f' 


(^2  +  ^2 

with  a  =  -(471/3), /i  =  -(p/ji),  of  are  odd  and  are  even  in  (p,  + 

€  iZ’®(R)  for  .s  >  0  with 


sup  (i|(l  +  (•,V’)||//>(R)  +  11(1  +  (•)^/’)IU^(R)) 

+  sup  (||^i^(-,t^)||//-«(R)  +  ||A^(-,l/>)||/f»(R.)) 

0<V^<oo  ^ 

+  sup  (||<^^i^(-,V’)IU'(R)  +  II<^A]'’(-,^)||/f.(R)) 

q<4!<Kx  ^  ' 

<  A'e, 


and  ii  is  a  constant  independent  of  e  but  may  depend  on  s  and  Ki. 

Since  if^(R)  is  embedded  in  bounded  continuous  function  space  (^^(R)  and 
satisfy  the  Cauchy- Riemann  equations  (20)  or  (24),  (1 -|-(^2y/2^±^  ^2 -f  (,o2)^/2a^  and  their 
derivatives  are  bounded  and  continuous  in  either  R  x  (  —  1, 0)  or  R  x  (0,  -f-00).  Then  from 
(34),  (35)  and  (40),  we  write 

f^iO  ,A  ,A‘^)  =  7ie0  +  , 


where  TZi  consists  of  all  nonlinear  terms  in  /2  and  p^TZi  G  iif*(R)  for  s  >  0.  Therefore 
(42)  can  be  written  as 


hikxtp)  F  [7ie(r7<p  +  ()] 


+  F 


h(k,rl^)F[n^)] 


(84) 


But  from  (45),  we  have 


'F[ir,^\  =  (-^^+pml  +  pAI=\h(k))-')  ’  (  biefl  + -F  [Kil 


Substituting  this  expression  into  (84)  and  taking  the  Fourier  and  inverse-Fourier  trans¬ 
forms,  we  can  obtain 


Up, 'll))  =F 


■k.jk,  ^)(l  +  7i(-7i  -H  p|A;|(l  -H  pe\kMk))-^)-^ )  F  [1l^] 

1-  ')iekz{k,  7/>)(l  -h  7i(-7i  p|A:|(l  -f  pe\k\k2{k))-'^)-^) 


However,  notice  that  71  <  0  and  k^{k^xl>)  is  bounded  and  even  in  k.  Since  'R.\  is  odd  and 
is  in  if^(R)  for  5  >  0,  integration  by  parts  twice  yields 


/'  +  =°  /  (fi  A  \ 

•'A)  =  sin(/j(^)  F  [T^i])  j  dk  . 

But  is  in  ff^(R)  and  the  derivatives  of  k4{k,ilj)  are  bounded,  which  implies  ip'^^  6 

ff'’(R)  for  ^  >  0  and  any  fixed  ?A-  Then  by  (45),  6  Therefore  from  (40) 

ip'^6~  €  H^(R)  for  s  >  0.  Now  we  use  integration  by  parts  in  (43)  and  the  oddness  of 
{l/ik)kzip{k,'tl})  to  have  (1  +  6  C'^(R)  for  integer  s  >  0.  From  (45)  and  the  fact  that 

^xp=o  ~  <9iU=o  =  0  in  (74),  ■ 


77  =F 


{  - 'Ti  +  +  P^lklhik)) 


-1 


A 

F 


p-\-CO 

J  ip 


77-2 


where  77-2  =  (3/2)(r;^)(^  +  /2,  and  F,2d(p  6  hr*(R)  using  Lemma  3.  Thus 

integration  by  parts  twice  yields  (1  +  G  (^^(R)  for  s  >  0.  Finally  we  use  the  Fourier 
transforms  of  exp(— |A;|e^)  and  e7(|A;|/fc)exp(— ejArl?/’)  in  (49)  and  (50)  and  the  classical 
estimates  of  elliptic  operators  in  [13]  to  obtain 

sup  ({|e+|  +  |A+|)(l+/  +  W))  <A' 

¥>eR.,V’€R+  '  ' 


and  the  derivatives  of  A"*"  are  bounded  by  K{1  +  .  Therefore  Theorem  5 

can  be  improved  as  follows. 

Theorem  6.  (32)  to  (35)  have  a  solution  stated  in  Theorem  5  with  6'^  and  A^  satisfying 
sup  {{\6f{p,i))\  +  \xf{p,tl^)\){l  +  p^  +  leV’l)) 

S 

+  J^sup  ((|Z)‘«f(v>,V’)l  +  ID'Af  (v,V.)|)(l  +  v"  +  k-AP))  <  A-£ 

J=1  ^ 

iox  any  integer  s  >  0,  where  is  denoted  as  the  f-th  derivative  with  respect  to  either  p 
or  xjj  or  both  and  K  is  independent  of  e  but  may  depend  on  s. 

After  transforming  back  to  9^,X^,  and  in  (30)  and  (31),  and  calculating  the 

Fourier  transforms  explicitly,  we  have 
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Corollary  1.  For  e  >  0  small,  the  equations  (20),  (21),  (24)  to  (26)  have  a  solution 


A'^(</’i,V’i)  = 

^~(<^i,’/’i)  = 


-2e'Vi«KeV’i  +/^)  ,2^+^  ,  ^ 


e2a^((ei/>i  +  A^)^  -  (e(/?i)^; 

((e,/., +/z)2  +  (e<^0^)- 

-2e-Vi«/W^(’/’i  +  1) 


(/x2+(e<^l)2)2 


+  (e(^i,’Ai), 


A  (</?!, •0l)  =  - 


ea// 


+  eAi  (e(^i,V’i), 


where  a  =  — (47i/3),yU  =  — (p/71),  and  Af  satisfy  the  properties  stated  in  Theorem  6. 

From  the  Bernoulli’s  equation  (6)  at  the  interface  ?/’!  =  0,  the  interface  of  two  fluids 
is  determined  by 


!,  =  (1  +  7,  e)-'(l/2)((l  -  )  -  p{l  -  e»^)) 


i>i=o 


= -A  +  eYi{e(pi)  =  eS{€(pi)  +  eY{e(pi) , 


where  sup^gjj^  |(1  +  (p'^)(d^Y (ip) / d(p^)\  <  Ke  and  K  is  independent  of  e.  Finally  we 

use  the  dimensional  variables  to  state  the  existence  of  the  solitary  waves. 

Theorem  7.  There  is  a  solitary  wave  moving  with  constant  speed 

U  ==  {{p~  -  p'^)9h/{p~{l  +  j^e))Y^^ 


in  a  two-fluid  flow  with  71  <  0  fixed  and  e  >  0  a  small  parameter,  where  the  height  of 
the  upper  fluid  is  infinite  while  the  depth  of  the  lower  fiuid  is  h,  and  the  densities  of  the 
upper  and  lower  fluids  are  p"*"  and  p“,  respectively.  The  profile  of  the  interface  between 
two  fluids  is  given  by 


y 


eau 


-2^3 


where  a  =  -(47i/3),p  =  -{p'^ / p  71), 


sup  1(1  -f-  x^){d^Y*{x)ldx^)\  <  Ke, 

and  K  is  independent  of  e. 


30 


Obviously  the  first-order  approximation  of  the  solution  given  in  Theorem  7  was  ob¬ 
tained  by  Benjamin  [4]  using  asymptotic  method. 
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